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This project explores the possibility of shallow ground water 
depth prediction from satellite data. Two methods are used for this purpose 
the computer compatible tape method and the Micro-densitometer methbd. ; 

Both methods are based on the fact that reflectance characteristics j 

■ ■ ■■■ ' . : . \ ' : f 

of a particular soil cover vary with variations in the ground water j 

'i, " . ' *’■ v’ *" * * /* g '* ' ■- .‘.i'’’’' , *'/ " 

depth. Calibration curves between reflectance units (CCT method) or ! 

light density (densitometer method) and water depth, are developed from j 

the field > ta and these curves are used for prediction of ground water 
depths. The CCT method is more accurate than the other method. It was found 
that these methods are not applicable for water depths over 10 m, or for 
hilly or built up regions. Very encouraging results were obtained in 



Kie prediction of ground water depths in alluvial soil C Gujrat 
region) in depths less than 10 m. For such regions this study 
establishes that the for macing ground water contour maps of large 
areas, the most economical and convenient method is a satellite 


survey® 
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CHAP 'TER - 1 


IN TR OPTI C TIGHT : 

1*1 Basic Principles ; The objective of this thesis is to explore 
the .possibility of predicting shallow water depths from satellite 
data. This study is aimed at providing a more economical and convenient 
medium (i.e. , a satellite survey) compared to ground surveys for 
making ground water maps of large areas. 

The area selected for study encloses a part of Rajasthan 
and a part of Gujrat (about 185 Km. x 185 Km.), A broad geological 
classification of this area is given in Pig. 1. The satellite data, 
pertaining to this specific area was obtained by LANDSAT in the form 
of a black and white photo negative , 

and Computer Compatible Tapes (CCT). This satellite data - the 
imagery and two CCTs - were procured by the Indian Space Research 
Organization (ISRO), from the EROS Space Centre, South Dakota, USA. 

Remote sensing is the science of obtaining k howl edge 
without physical contact. The remote sensing satellite by the use 
of Multi-Spectral Systems (MSS) can gather such 'knowledge' from an 
altitude of several hundred kilometres. (920 Km. for LANDSAT), by 
scanning the land masses underneath and recording the neglected 
solar radiation (or 'reflectance') in different wave-bands. For 
water studies, we use. the reflectances recorded by the MSS in the 
infra- red region or band - 7 (0. 8-1.1 m. wave- length) because In 
this band water appears black on the positive print due to strong 
absorption bands in the near Infra-red part of the spectrum. Eor a 
particular soil in a given area, the variations in shallow water 
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depths alter the tone of the ground surface and these variations, 
most discernible in Band - .7, alter the reflectances of the same 
ground surface depending upon the proximity of the ground water 
to the ground surface. Hence che variation in reflectances, • 

of the same ground surface is related to the variation in the ground 
water depths. By measuring the ground water depths by bore holes etc. 
at a number of distributed points in the given area and by recording 
the corresponding reflectances from the available satellite data, 
we can study the correlation and plot a curve of reflectance versus 
ground water depth. This is a calibration curve, and it can also be 
called a calibration - prediction curve because it can be used to 
predict the ground water depths at different points in the given 
area knowing the corresponding reflectances. 

Two techniques based on this principle have been employed 
in this research study for the prediction of shallow ground water 
depths. One uses a. micro-densitometer and the other is' based on the 
application of magnetic tapes (CCTs). The COT approach required 
an extensive use of the DEC-10 computer available in IIT Kanpur. 

The promising results obtained in this study highlights 
the scope of satellite application for predicting shallow water depths. 
For a large area survey like a state survey or a national survey 
for marking ground water depth contours, it is obvious that the 
most economical and convenient option is a satellite survey. 

"T* -4* 






CHAPTER - 2 
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MAGNETIC TAPE METHOD ; . ' 

2.1 Computer Compatible Tape (CCT) Approach 

The data collected from a Multi spectral scanner (MSS) 
on a aircraft or satellite can he stored on a Computer Compatible 
Tape after being processed at a ground data collection centre. The 
MSS system collects ground reflectance data by scanning along lines 
perpendicular to its flight path and recording reflectance values 
as 'pixels’ along each 'line'. The resolution of a pixel in La.ndsat-3 
is about 80 m. After being processed at the ground station to 
eliminate undesirable interferences or 'noise', the data or 'signals' 
is recorded on a CCT in the form of ’records'. The 'pixels’ or reflectances 
are recorded in 'bits' on a. CCT in the binary mode. 

2.2 CCT Storage Mode s 

Information is stored in a CCT in the form of 'records'. 

Each record represents one line scanned by a satellite in a given 
imagery. Each record is subdivided into 'words' - a smaller storage 
unit, i.e. And each word is subdivided into 36 'bits', a bit being 
the smallest storage space on a CCT. Each bit stores one binary 
number. In our case, one word ' in the CCT Stores four pixels, each 

pixel occupying seven bits of space (hence the maximum reflectance 

7 

value cannot exceed 2 = 128). The reflectance values range from 
0 to 128. 

As the MSS system records data in Bands 4, 5, 6 and 7. 

The same pixel is recorded on a '800' density or'l600 # BPI tapes 
in four different parts in the four wave -ands. Each '800' density 
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CCI contains two sections or 'files', and each file contains data 
of the sane one imagery in a given wave hand. Thus to record the same 
imagery in Bands 4, 5 , 6 , 7 ately, two '800' density CCTs are 

required, each -wave hand confined to one file on a. tape. 

The CCT containing Band - 7 part of the imagery was used 
for ground water study in this study because variations in ground water 
profile are most discernible in Band - 7 (0. 8-1.1 pm wavelength) j 

A Computer program (see Appendix No. I ) was developed 
which gives the line number and the pixel number of a particular spot 
on the imagery or a. band map of the given Rajasthan - Gujrat region 
for the given latitude and longitude coordinates of this particular 
spot.. Hence, to read the reflectance value 

from CCT at any given point on a imagery, first the latitude and 
longitude of this point are determined from the topo sheet and these 
coordinates are fed into the computer program to get the required 
line number and pixel number of the given point. Once the line number 
and pixel number are known, the CCT is used to read the records on the 
tape till we reach the required record (or line scanned by the MSS) 

Next we activate Computer program (see Appendix II) which reads 

the record, converts the data from binary to decimal mode and prints the 

given reflectance (or 'pixel' scanned by the MSS). 

2 • 3 C en t ral pixel Meth o d & Grid Averaging Method s 

Four sets of correlation were attempted between reflectance 

the so 

and ground water depth. In the first esse or / .called the 'central 
pixel method', only the exact pixel, corresponding to the given latitude 
and longitude coordinates of a given point on a. land map or imagery, 
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was oelec.ed and its reflectance value was correlated with the 
corresponding water depth. 

In the arid Averaging Method, we take a number of pixels 
on either side of the central pixel and take the arithmetic mean of 
the reflectances of these pixels which form a square or grid with the 
central pixel occupying the centre position. It is obvious that for 
U oquorc gnu, each side/the grid must have an odd number of pixels. 
Hence we can have a grid with 3 pixels on each side, 5 pixels on each 
7 pixels on each side etc. In this study these three grid values 
cnoocn (3, 5, 7 pixels, i.e.). T t can easily be seen that in the 

caoo of the grid approximation method for 3, 5, 7 pixels making up the 

id sides we will hpv*p +• 4.ni rr , 4_ ■u . 

nave to take the arithmetic mean of 9, 25 and 4-9 

reflectance valves respectively to get the representative reflectance 

gno. s' or a large number of control pixel points, manual 

computations being time consum ing for grid averaging at each point, 

3 S “ Ple 00m P uter P^am was used which gave the representative 

reflectances tor both central pixel method and grid Averaging Method 

(o. 6, 7 pixels grids) far each observation point. 

* 

After this Step is completed, knowing the latitude and 

longitude of a given observation point, four representative reflectances 

corresponding to this observation point (one for the central pixel , the 

other three for the 5, 6 and 7 pixels sise grids) are computed. The 

next step is the preparation of rcfloct-inrd , 

reflectance versus depth •'•curves* ' • 

2.4 Calihrati on- 

The four representative reflectances, obtained from the 
previous step, were plotted individually against observed ground water 
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correlations were attempted. After the points were marked on each graph 
the next step was curve fitting. A qualitative inference about the 
nature of the calibration curve fojfc the magnetic tape approach,- was 
deduced from basic principles. As water appears black in Band ■ 7 the 
shallower the 'water table is the more affected the ground surface 
would be from the moisture proximity. 


Hence for lesser water depth the same ground surfa.ce will 
give lesser reflectance and for deeper depth the reflectance will be 
greater. For a surface water body like a lake, the reflectance will be zero, 
the minimum reflectance value. Hence we have an exponential growth 
curve of increasing slope, between reflectance values and the corres- 
ponding water depths. Mathematically, the equation of such a calibration 
curve fitted is s- 
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bx+c 

e 


(II) 

where 1 b’ and ' c’ are constants, 
y = reflectance, 
x = ground water depth 


In order to fit a calibration curve obeying Equation (II) 
in the four sets of graphs between the representative reflectances and 
the observed ground water depths, the method of least squares was 
ad opted . 


Evaluation of t b' and ' c \ in - Equation ( II ) ; 

Taking the natural logarithm of both sides in Eqation (II), 
log e y - hxfc 

let Y = log e y, b = BB, c = cc 
Y = BB X + CO 



It can easily be established that the normal equations are, 

= 33 Tx + n CC - - - (III) 

X.xY = BB^x' 2 + CC& - - - (IV) 

where n = number of observed points. 

The values of ’ BB’ and ’CC’ can be found by simultaneously 
solving Equation (III) and (IV) b^ r inserting the observed values of 'X' and 
' Y ’ . 

Once ’ BB' and ' CC ' are known the standard error, , for the "fit" 
determined froms- 



After the 3-S* test 
correlation was determined to 
between ’x' and T y * and if so, 
coefficient is defined as: 


had been applied, the coefficient of 
find whether there was any correlation 
how good it was. The correlation 

J(x- x) (y - y) 

■yf J (x-x ) 2 £ (y-y ) 2 


The value of 'r ' can vary from -1 to +1. A negative ’r ’ indicates 
that ’y* decreases with increase in 'x*. If r - 0, it indicates an absence 
of correlation. If r 0.5, a weak correlation exists? if r-i.0.8 - 1.0, 
a. good correlation is concluded (subject to the usual standard statistical 
tests with the given level of confidence). 

Results s In this study, with the help of a computer program (see 
Appendixm), the calibration curve equations, standard deviations and 
coefficients of correlation were determined for the four plots in the 
magnetic tape method for the Gujrat region. These results are given in 
Table - 1 . 
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Table - 1 


! “ 1 

? 

! 

? 

? 

} METHOD 

? 

V 

- • 

— j r— 1 ~ t 

CALIBRATION CURVE EQUATION ! STANDARD \ CORRELATION j 

BBx+cc ’ DEVIATIONS; COEEEICIENT * 

V = c ; i , : 

BB 

CC J f or 1 x ' 

! 

for 'y| ’r* J 

r 

t 

? 

i 

jCentral Pixel 

1 

? 

*3 pixels grid 

? 

f 

l5 pixels grid 
? 
f 
f 

'7 pixels grid 

f 

t 

0.20175720 

0.16010590 

0.07443128 

0.03375640 

— — 

! 

t 

1.92411720 ; 2.047 

1 

2.41294530 1 2.047 

1 

3.08453320 j 2.047 

t 

3.46898590 < 2.047 

! 

1 

t 

" T“ — f 

; j 

0.44& 0.9389 j 

! I 

T \ 

0.36^0.8916 j 

I ; 

0.179J 0.8523 j 

J 1 

t ? 

I 0.125} 0.5544 j 

t ! ? 

t 1 , 

t > . t 


The field data for the Gujrat region used for obtaining these, 
results are given in Table 2( Appendix IV). 

The field data for the Rajasthan region - used in making 
reflectance vs. depth plots (see g ig."' - ®v^ ; by central pixel method and 
3 pixels/ and 5 pixels linear averaging methods - are given in Table 3. 

No appreciable correlation was found between reflectance 
and depth in the Rajasthan region after inspecting figs. 2-3. 


Table - 2 


Field data Gujrat region (CCT Method ) 


sia r 

no., 1 Longitude! Latitude } Line 
{(Degrees ) { (Degrees ) {No, 

» i i 


p , ,J Reflectance 

No' Centra ^ 3 -pixel} 5 -pix.J 7 -pixel } Depth 
°* [ pixel { grid j grid 
5 method f ! 


1, 

72,850000 

24,340000 

1436 

2195 

19 

27 

31 

35 

5,0 

2, 

72,425000 

23.915000 

2163 

1764 

59 

57 

45 

42 

10,0 

3. 

72,366666 

24.266666 

J656 

1524 

21 

24 

31 

35 

5,8 

4. 

72,366666 

23.916666 

2174 

1679 

34 

43 

42 

48 

* 

CD 

5 . 

72.020000 

24,520000 

1357 

910 

27 

26 

32 

39 

6,8 

6. 

71,810780 

25.354423 

172 

241 

58 

59 

47 

43 

9.9 

7. 

71.866667 

24.300000 

1721 

781 

24 

39 

41 

45 

7.7 

8. 

71,625000 

25.363000 

437 

167 

19 

22 

30 

36 

4.8 

9. 

71. < « •?!%«* 

24,06 1374 

2048 

1171 

24 

37 

40 

47 

5,4 




Table — 3. 
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Field data .*»' ■■ Ra j a.sbhan region : 

( Longitudes and Latitudes of these field points are given 
in Appendix V.) 


s5“ 

/ No* 

1 " 

| Line 

l NO. 

5 

1 Pixel 
| No* 

1 7. HEFLECTANCF 



J ' “ • 

5 Central I 
| pixel l 

3 pixels 
average 

I 5 pixels 
| average 

5 7 pixels 

5 average 

* Measured 
! depth 

I (m.) 

It 

179 

242 

63 

47 

44 

42 

9,90 

2f 

235 

920 

22 

41 

44 

47 

3,46 

3t 

236 

1320 

66 

59 

50 

49 

8,64 

4. 

248 

1256 

4i 

42 

43 

49 

6.23 

5. 

261 

1179 

25 

44 

47 

48 

8.66 

6 # 

271 

1296 

50 

51 

50 

42 

2.11 

7. 

278 

764 

34 

30 

41 

41 

5.72 

0. 

205 

12 00 

59 

45 

42 

41 

8#45 

9. 

206 

025 

27 

40 

45 

47 

6,90 

10. 

299 

1200 

39 

39 

a? 

44 

0.16 

11. 

320 

607 

62 

52 

50 

47 

5,90 

12. 

337 

609 

35 

31 

42 

46 

2,99 

13. 

330 

776 

30 

39 

39 

45 

5.69 

14. 

340 

644 

40 

52 

43 

39 

7,72 

15. 

351 

049 

60 

63 

51 

46 

3, SO 

16. 

305 

935 

90 

29 

36 

39 

6,27 

17. 

419 

630 

45 

41 

SO 

32 

3,00 

10, 

426 

1136 

63 

57 

55 

47 

6,30 

19, 

' 436 

003 

22 

33 

36 

42 

8.27 

20. 

430 

167 

27 

31 

44 

47 

4,80 

21, 

447 

185 

61 

62 

50 

4$ 

6,39- 
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Contd, Table — 3 


si. { I : j " ' 'REFLECTANCE " 1 “ 7 la 

No. { Line I Pixel J Central I 3 pixels {5 pixels p7 pixels I deoth^ 

■ l No; l No, i pixel { averavge I average {average J f \ 

l L__l I _i I ! ^ 


22. 

455 

.1035 

6-9 . 

. .. 60 ■ ' 

58 

53- . 

0*75 

23. 

456 

754 

56 

46 

42 

• 40 ' 

2*77 

24. 

460 

675 

54 

49 

40 

41 

3 # C6 

25. 

■467 

1203 

' '63 

- • 60 

50 

52 

0*90 

26, 

403 

909 . 

26 

43 

52 

55 

7.00 

27, 

409 

73 4... 

65 

61 

49 

46 

2*00 

20, 

479 

631 

62 

64 

53 

44 

4.75 

29. 

5 09 

1237 

31 

32 

33 

39 

6.45 

3 3, 

519 

1177 

19 

29 . 

33 

34 

9.2 0 

31, 

525 

724 

25 

3 0 

42 

49 

9.95 

3 V 

533 

1160 * 

30 

36 

36 

30 

6.5 C 

33* 

553 

737 J 

41 

36 

51 

56 

3*67 

34* 

563 

1232 . 

53 

47 

47 

40 

5.09 

35, 

629 

435 

• . 24 

36 

42 

46 

6,10 

36 . 

646 

692 

27 

30 

43 

40 

5 *1 C 

37. 

3 09 

216 ” 

54 

41 

30 

37 

3.2 0 

on 

Ou* 

903 

414 

50 

57 

45 

42 

9,53 

*9. 

1251 

314 

9 

15 

20 

23 

4*90 

4 • j * 

136 0 

9 09 

23 

35 ’ ' 

30 

39 

6.0 0 

41. 

1602 

4 02 

62 

53 

40 

49 

4.0 0 

42. 

2 C 47 

776 

34 

37 

37 

43 

7.90 

43. 

2127 

4S4 

35 

32 

44 

50 

0.40 



Reflectance v». Depth plot 

Ji tsX ♦ CPX 

( R*j Mthan region ) 






Reflectance v*. Depth plot ? : 'jV j 

£**•-£ CCT 

-. i . ; ; (Rajasthan rafian) i 







Reflectance vs. Depth plot 



CCT method 


L££L*!*« ) 

* : (Rajasthan region) 



4.0 $.0 6.0 • - v ' 7.0 . ' 8,0 . ' 8.0 • 10 0 

^ ^ •' ■ 
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2. 5 Prediction Process (Magnetic Tape Method ) ; 

Prediction Method A number of points were taken on a 1” = 4 miles 
scale ground water contour map of the given region and their latitudes 
and longitudes were determined. Using these latitudes and longitudes in 
the ' 1ATPIX' computer program (see Appendix I), the line numbers and pixel 
number corresponding to each point were found. Then using these line 
numbers and pixel numbers and the 'TAPE. FOR' computer program (see ^ • • • 

Appendix II) the corresponding reflectances (pixel values) were found. 

Next with the help of the . calibration curves and these reflectances 
the corresponding ground water depths were predicted for the given set 
of points taken on the man. These predicted ground water depths were 
then compared with those indicated in the conventional ground water 
contour map to complete the prediction process. 

To attempt a prediction of ground water depths in the G-ujrat 
region ( lower half of the given imagery) the Banaskantha district which 
mostly" covers this part was selected. This district area was subdivided 
into three sub regions (see Pig. 15)for analysis. Important details about 
each region are given below; 

Region (l) Mostly alluvial with ground water depths < 10 m 

LcJc » “t 

longitude 71.45 - 71.85, la-tidue s 24.00 -,24.50 

Region (2) Alluvial transition to hilly region, water depths ,p 10 m 

longitude i 71.90 - 72.30^ latitude : 24.05 - 24.50 

Region (3) Mostly hilly region with alluvial regions in minority. 

Depths s variable? includes > 1 Op, depths also. 

longitude i 72.35 - 72.70'', latitudes 24. 10’- 24.25' 

Region (l) was expected to be most suitable for prediction 
purposes and the results of the two prediction methods - tap.e method & densitc 
meter method bore this out. In Region (2) & Region (3) the methods are 
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unreliable because of the deep depths mid hilly terrain# 

MM£J1 Ml JMMJMSM fh * *l©ld data etc# isgiven in table* 4-7* 

A separate case study was taken up to show that the satellite 
survey methods are not useful for detecting depths deeper than 1C a# 
Uinetton field measurements of water depths greater than 10 m were 
undertaken and the predict i m nroodss using the tap© method was 
applied { 8©e table — ? )* ^e following result® were obtained s 

H m 19 points# O' « 2 #75 a 

Hie 3 - O’ test sheas that we have to reject any x> 8*27 a 
hence we roject serial no# 16 obaerra* ion taHe - 7* '"herefor-e, we 
finally get t 

ifcia, 

Standard error, CT * 2*63 ® 

flic high value of the standard error proved the point 
adequately that satellite survey methods are unrelldtle in depth© 
greater than 10 »# 

team, m * ( See table - 4 ) 

Humber of observation points* I 61 
Standard error* (f n • t, 1 # 36 a 

The 3 - O' test ©how© that we have to reject any x ?- 4*08 a 
hence we reboot the following two observations » x * 7#84 a {Longitude 
71 * 6 °* latitude * 24*4° ) and % =r 6.13 a C Longitude m 71 # 75 °* 

Latitude « 24.15° )• 

After removing these two point© we get I 

1=59* <r n rn + 0.99 a* 
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The standard error is t 0*99 a* The tape method thus fives 
a food result in He ion (1), which ismostly alluvial In nature and 
with water depths less than 10 m t where our calibration - prediction, 
equations yield the best results# 

Horion (2) s { See table ** 5 5 

N « 61 (T* 4*23 a 

'The 3 *• <T test shows that v;e have to reboot a ry x>12»69 is* 
hence we rejoot the tea observa-* ions corresponding to serial nos 28* 
31* 36, 37, 36* 42* 43* 46* 48 and 52 in table * 5* 

After r amoving these ten point® we get * 

K « 51* 

Standard error* 0" * 3*39 a 

Therefor©* the standard error being 3*39 m# we can conclude 
that in Region (2) we get highly unreliable results# This region 
is mostly of more than 10 m ground water depth* where our satellite 
survey methods ore not applicable# 

IMsaJlLi < Sfi « taw* - 6 5 

s *■. to* <T* 1*97 m 

The 5 • (Ttest shows that m have to reject any x > 5*91 m* 
so we reject one observation correspondin'; to eeridl no* 10 in tttt 
table - 6* 

therefore, we have * 

JST ** 9* 

Standard error* <f * 1 *84 m 

Y/e observe that a few points ( for ©at** serial no* 3 la 
table » 6) show large differ eases between predicted ant interpolated 
depths* 
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Mb occurs most probably because the nature of the local 
terrain, which is hilly (extension ©f Mount Aha hill region)* A few 
points show a tea sgr good agreonenfe between interpolated end predicted 
depths (for ex* # serial no* 9) and this could he attributed to allurial 
patches in between the part of the hilly regions for ground water depths 
less than 10 m* 



Table — 4 
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Field data & calculations for prediction 
(CCT method! Region (1 ), Gujrat) 


SI. 

No. 

i : i 

j Longitude | 
l (Degrees )1 

Latitude 

•(Degrees) 

I Inter- 
| polated 

I Depth, D 
l (m. ) 

1 

J Reflect— 

{ ance 

I 

"1 

i Predicted 

v depth, P 

I (m.) 

5 

I 

5 

Deviation 
x = (d-p; 
(m. ) 

1. 

71.45 

24.20 

5.00 


27 

6.80 


1.80 

2. 

71.45 

24.05 

9.00 


25 

6.42 


2.58 

3. 

71.45 

24.00 

9.00 


59 

10.67 


1.67 

4. 

71.50 

24.25 

5,00 


22 

5.78 


0.70 

5. 

71.50 

24.20 

6*40 


33 

7.79 


1.31 

6. 

71.50 

24.10 

6.81 


36 

0.23 


1.42 

7. 

71.55 

24.30 

4.00 


23 

6 * 00 


2.00 

0. 

71.55 

24.20 

6.67 


58 

10.59 


3.91 

9. 

71.55 

24.15 

7.00 


50 

10.59 


3.59 

10. 

71.55 

24.10 

7.00 


51 

9.95 


2.95 

11 . 

71. 5S 

24.05 

3,59 


61 

10.84 


2.25 

12. 

71.55, 

24.00 

8.10 


32 

7.64 


0.54 

13. 

71.60 

24.40 

3.00 


61 

10.04 


7.84 

14. 

71.60 

24.35 

4.33 


24 

6*22 


1.04 

15. 

71.60 

24.30 

7.00 


2® 

6.00 


0.20 

16. 

71.60 

24.25 

9.00 


72 

11.66 


2.66 

17. 

71.60 

24.20 

7.00 


24 

6.22 


0.88 

13. 

71.60 

24.15 

7.00 


26 

6.61 


0.39 

13. 

71,60 

24.05 

7.40 


30 

0.49 


1.01 

20. 

71.65 

24.05 

7.06 


26 

6.61 


0.45 

21. 

71.65 

24.45 

5.64 


45 

9.33 


3.69 

22. 

71.65 

24.40 

4*8 0 


25 

6.42 


1.62 

23. 

71.65 

24.25 

7.00 


22 

5.78 


1.22 

24. 

71.65 

24.20 

7.00 


40 

9.65 


2.65 

25. 

71.65 

24.15 

7.00 


33 

7.77 


0.77 

26. 

71.65 

24.10 

7.00 


30 

8.49 


1.49 

27. 

71.65 

24.00 

5.03 


20 

5.70 


0. 05 
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Contd. Table — 4 


Si . 

No. 

{Longitude 1 Latitude I 
$ (Degrees ) | (Degrees ) | 

I l l 

Inter— I 

polated { Reflect— 
Depth, D 5 ance 
Cm .) I 

f 1 ? re S ct ® d ! Deviation 

0 depth, P i 

1 (m.:) ' ix = 5D-P! 

5 ' I (m.) 

20 • 

71.70 

24.50 

7.44 

27 

6.60 


0.74 

29 . 

71.70 . 

24,45 

7.29 

20 

5.31 


1.98 

3 0 . 

71.7 0 

24.40 

6.47 

25 

6.42 


0 . 05 

31 . 

71.70 

24.25 

7.00 

52 

10.05 


3.05 

32 . 

71.70 

24.20 

6.00 

20 

5.31 


0.77 

33 . 

71.75 

24.50 

0.02 

21 

5.55 


2.47 

34 . 

71.75 

24.45 

, 7.77 

29 

7.15 


0.62 

35 „ 

71.75 

24.40 

. 7.51 

52 

10.05 


2.54 

36 . 

' 71.75 

24.35 

7.24 

17 

4.51 


2.73 

37 . 

" 71.75 

24.30 

7.00 

45 

9.33 


2.33 

33 . 

71.75 

24.25 

7.00 

27 

6.00 


o^Z 0 

39 . 

71.75 

24.20 

. 5.00 

26 

6 . 6 * 


1.61 

40 . 

71.75 

24.15 

. 5.17 

■ 67 

11.3 0 


6.13 

41 . 

71.75 

24*10 

, 5.64 

25 

6.42 


1.34 

42 . 

71.75 

24.00 

5.00 

29 

7.15 


2,15 

43 . 

71,00 

24.50 

0.65 

25 

6.42 


2.25 

44 . 

71.0 0 

24*45 

0.40 

44 

9,22 


0.74 

45 . 

71.00 

24.40 

0.20 

32 

7.64 


0.64 

46 . 

71.80 

24.35 

0.00 

25 

6.42 


1.50 

47 , 

71,00 

24,30 

9.25 

37 

0,36 


0.09 

4 G . 

71.00 

24.25 

9.63 

30 

8.49 


0.14 

49 . 

71.00 

24,20 

9.00 

49 

9.75 


0.75 

5 0* 

71.00 

24.15 

7.00 

40 

0.75 


0.95 

51 . 

71.00 

24,05 

7.40 

30 

0.49 


1.01 

52 . 

71.00 

24.00 

5.10 

33 

7.79 


2.60 

53 . 

71,05 

24.05 

9.24 

50 

9,05 


0.61 

54 . 

71.05 

24.45 

9.20 

26 

6.61 


2.67 

55 . 

71.05 

24.40 

9.16 

20 

6.90 


2.18 






Table — 5 


Field data & calculations for prediction 
( CCT method l Region (2)/ Gujrat ) 


SI. 

NO. 

1 

i Longitude 

j (Degrees) 

$ 

Latitude 

5 (Degrees) 

I 

| Inter - * 

« polated x 
* Depth, D I 

1 (m.) | 

Reflect- 

ance 

I Predicted 

8 depth, P 
| Cm, ) 

1 

j Deviation 
$ X = ;D — P { 

i (m.) 

■1. 

71.90 

24.50 

9,90 

24 

6,22 

3.60 

2. 

71.90 

24.45 

10.14 

23 

6.00 

4.14 

3. 

71.90 

24.40 

10.50 

53 

10.14 

0.44 

4. 

71.90 

24.35 

11.77 

32 

7.64 

4.13 

5., 

71.90 

24.30 

13.75 

30 

0.49 

5.26 

6. 

71.90 

24.25 

14.76 

25 

7.42 

0.34 

7. 

71.90 

24.20 

14.76 

62 

10.92 

3.04 

0. 

71.90 

24.15 

14.00 

76 

11.93 

2.07 

9. 

71.90 

24.10 

15.0 

32 

7.64 

7.36 

10. 

71.90 

24.05 

11/40 

61 

10.04 

. 0.56 

11. 

71.90 

24.00 

7.50 

30 

0.49 

0.99 

12. 

71.95 

24.50 

10.53 

29 

7.15 

3.30 
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Contd. 


Table 


31. T Longitude 1 Latitude 1 
No. $ (Degrees) 5 (Degrees) Jj 

I 5 I 

Inter — jl _ ' , 

polated ? Ref lect “- 
hepth, D | ance 
(m* ) 1 

1 Predicted 

5 depth* P 
j ( m.) 

I Deviation 

x = ' D—5 1 . 

1 (m.) 

13. 

71.95 

24,45 

10,75 

47 

9*55 

1.20 

14. 

71.95 

24.40 

12.33 

31 

7.40 

4.05 

15. 

71,95 

24.35 

12.00 

30 

0*49 

3,59 

16. 

71.95 

24.30 

14*27 

10 

4*79 

9.40 

17. 

71,95 

24,25 

19*33 

65 

11.15 

6.10 

10. 

71,95 

24,20 

17.33 

22 

5.70 

9.10 

19. 

71.95 

24.15 

17*27 

36 

0.23 

9.04 

2 0. 

71.95 

24,10 

14.50 

59 

10,67 

3.05 

21. 

71.95 

24.05 

10.25 

27 

6.00 

3.45 

22, 

71,95 

24.00 

10.20 

62 

10.92 

C.72 

23. 

72.00 

24.50 

15.00 

47 

9.55 

5.45 

24. 

72.00 

24.45 

13.00 

26 

6.61 

6.39 

25. 

72.00 

24.40 

12.50 

42 

0.99 

3.51 

26. 

72.00 

24.35 

14.20 

17 

4.51 

9.69 

27. 

72.00 

24.30 

10.57 

27 

6.00 

11.07 

20. 

72.00 

24.25 

19.70 

14 

16.16 

12.03 

20. 

72.00 

24.20 

19.66 

36 

0.23 

11.43 

O 

W m 

72.00 

24.15 

17*67 

64 

11.10 

6.56 

31. 

72.00 

24.10 

14,40 

0 

0,77 

13,63 

332. 

72.00 

24.05 

10*50 

45 

9.33 

1.17 

33. 

72.05 

24.50 

19.00 

29 

7.15 

1.1,05 

34. 

72.05 

24.45 

19.00 

27 

6., 00 

12.30 

35. 

72.05 

24.35 

21.00 

46 

9.44 

10,66 

36. 

72.05 

24.30 

22.00 

26 

6,61 

15.30 

37. 

72.05 

24.25 

22.00 

26 

6.61 

10.39 

33. 

72.05 

24.01 

17.00 

13 

3.10 

13.02 

39. 

72.05 

24.10 

14.27 

60 

11.30 

2.89 

40. 

72.10 

24.40 

19.00 

• 26 

6,61 

12.39 

41. 

72.10 

24,35 

19.33 

66 

1 0. 42 

0.91 

42. 

72.10 

24.30 

20.67 

27 

6.00 

13,97 

42 , 

72.10 

24*25 

22.61 

36 

0.23 

14,39 


Contd 


Table — i 
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Si. 

No. 

7™ “ 

5 Longitude . 
J (Degrees ) 

5 

1 

' $ Latitude 

5 (Decrees) 

5 

1 Inter — 

5 pointed 

1 Depth, D 

1 (m. ) 

T 

5 

5 

5 

Reflect - 

ance 

5 Predicted 
l depth, P 
i (m.) 

| 

5 Deviation 
lx = ‘D -Pi 

5 ‘Cm. ) 

JS_ 

44. 

72.10 

24.15 

16,70 


76 

11.93 

jT 

4.06 

45. 

72.10 

24.10 

14.67 


26 

6,61 


0.06 

46. 

72.15 

24.35 

19.00 


23 

6,00 


13* 00 

47 * 

72.15 

24.30 

19,40 


49 

9,75 


9.65 

40. 

72.15 

24.25 

19,90 


22 

5.70 


14.12 

49. 

72.15 

24.20 

10.75 


26 

6.61 


12.14 

50. 

72.15 

24.15 

16,05 


50 

.9.05 


6.19 

51. 

72.15 

24.10 

14.53 


36 

0.23 


6.31 

52. 

72.2© 

24,3 0 

19,00 


19 

5.06 


13.94 

53. 

72.20 

24.25 

19.00 


53 

10.14 


0.06 

54. 

72.20 

24,00 

10,33 


66 

11.23 


7.10 

55. 

72.20 

24.15 

16,14 


35 

0.09 


6.29 

56, 

72.20 

24.10 

JML.ll 


66 

11.23 


2.00 

57. 

72.25 

24.25 

17.50 


53 

10.14 


7.36 

50. 

72,25 

24.20 

10.23 


34 

7.94 


10.29 

59. 

72,25 

24.15 

15.95 


47 

0.36 


7.59 

60. 

72,30 

24.25 

13.00 


23 

6.00 


7.00 

61. 

72.30 

24.20 

16,11 


59 

10.67 


5.44 
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Table — 6 


Field data & calculations for prediction 
( CCT method * Region (3), Gujrat ) 


si. 

NO. 

T“ 

r Longitude 
| (Degrees ) 

1 

-5 — - 

l Latitude 
| (Degrees) 

| Inter — J Reflect 
» polated I ance 
j Depth, D I 

1 (m. ) I 

— 1 Predicted 

5 depth, P 

5 (m.) 

5 

1 Deviation 
l X =-D -P : ' 

I tiu ) 

1. 

72.35 

24.25 

7.57 

44 

9.22 

1.65 

2. 

72.35 

24.20 

x 0.17 

20 

5.31 

4.06 

3. 

72.40 

24.25 

5.00 

59 

10,67 

5.67 

4. 

72.40 

24.20 

5.00 

29 

7.15 

2.15 

5. 

72.45 

24.15 

9.80 

26 

6.61 

3.19 

6. 

72.50 

24.15 

9.00 

50 

9.05 

0.05 

7. 

72.55 

24.15 

11.00 

29 

7.15 

3.05 

8. 

72.60 

24,15 

11.00 

20 

6.90 

4.20 

9. 

72.65 

24.15 

11.00 

59 

10.67 

0.33 

10. 

72.70 

24.15 

11,00 

19 

5.06 

5.94 
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Table - 7 

Field data & calculations ( CCT method .Gujrat region; 




for 

depths >10 m. 

) 




SI,. 

No. 

* Latitude 
| (ijegrees) 

T" 1 

x Longitude 
* {Decrees) 

Measured 

I depth 

5 (m. )D 

Line 

l HO - 

1 Pixel 
j NO. 

Reflect 
* ance 

1 

5 Predicted 
$ Depth (m)P 

l 

Deviation 
$ x= * P— D ; 

1 Cm. ) 

■i; 

24,280. 

.< 72,180' : 

11.5 0 

1679 

1247 

51 

9.95 

1.55 

2. 

24,220' 

72,020 " • 

12.20 

1804 

1040 

48 

9.65 

2.55 

3. 

24.350 

72.530 "0 

10.80 

1495 

1725 

51 

9.95 

0.85 

4 , 

24.190 , 

72.430' 

11.20 

1755 

165 0 

59 

10.59 

0,61 

5. 

23,875 . 

72.635' 0 

12.00 

2173 

2089 

48 

9.65 

2.35 

6 . 

23.850' 

72,640' • 

12.00 

22 09 

2107 

27 

6.70 

5.30 

7, 

23.825 

72.645' ' 

12.00 

2245 

2125 

17 

4.51 

5.49 

8, 

23.85 0< 

72.665 

14,00 

2203 

2143 

36 

8.23 

5.77 

9. 

23,825 

72.675' " 

14.00 

2238 

2169 

24 

6.02 

7.78 

10 , 

23,800' 

72.680. 

14.00 

2274 

2188 

45 

9.33 

4.67 

11, 

23,85 0 

72.690' * 1 

16.00 

2197 

2180 

33 

7.79 

8.21 

12 

23.825 

72.698. ’ ' 

16.00 

2232 

2203 

33 

7.79 

8.21 

13, 

23.800' 

72.700 

16,00 

2269 

2217 

50 

9.85 

6.15 

14, 

23.825 r 

' 72.727' ' 

18,00 

2226 

2245 

55 

10.33 

7.67 

15, 

23.800' " 

: 72,727 

18,00 

2262 

2257 

52 

10.05 

7.95 

16, 

23.772' " 

' 72.737 

2 0.00 

23 01 

2284 

59 

10.67 

9.33 

17, 

'54,266 

72,183 : 

14.50 

1698 

1257 

45 

9.33 

5.17 

18- 

24.174 1 . 

- 7 2,450 

10.80 

1770 

1698 

35 

8,09 

2.71 

19. 

24.081' 

. 71,941 ' . 

10.00 

2 027 

986 

26 

6.61 

3.39 
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CHA? 'TER - ^ 

MICRO DEN S I TOMB TER METHOD 

3 • 1 Micro Densitometer Approach i 
3.1.1. General s 

The micro-densitometer is an instrument which measures the 
light density at any point on the negative of a photo imagery. The 
Transmittivity , T, is defined as the ratio of the amount of light 
transmitted to the negative to the amount of light incident on the 
negative at any particular point on the negative. So ’ i" can vary from 
’ 1 1 (transparent part of the negative) to ’O’ (Opaque part of the 
negative), light Density, D = log 1Q (— ~) (1). 

Variations in surface features on land are reflected as 
variations in the shades or tones of a Black and White negative of an 
imagery. These variations in shades can he expressed in terms of 
variations in light density. Hence light density is an index of the 
properties of surface features mapped in an imagery. The microdensitometer 
gives o. plot of transmittivity, T, versus distance along any straight 
line etc. on a photo - negative. The value of ’ T 1 at any desired spot 
can be noted from this plot and the corresponding light density -'D* can he 
found cut from equation (l). 

3.1.2 Case Study : 

Microdensitometer plots for this imagery were obtained from 
the Satellite Applications Centre, Ahmedabad (Gujrat). The photo negative 
was divided into six equi spaced lines such that three lines covered 
the Rajasthan region the other three covered the Gujrat region. 
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Continuous Micro Densitometer plots for each of these lines were then 
obtained . 

3 . .1 . 3 Calibration :Part : 

A field study was undertaken to the Rajasthan and Gujrat 
regions and 50 wells were mapped and there water depths the latitudes 
and longitudes of the nearest land mark on the map noted, a.s well as the 
orientation and the distance to the landmark from each well. The corrected 
latitudes and longitudes of each well were then determined from the 
following relation: 


«(. = Angle between the north direction 
and the.|^g to the Actual location 
from the assumed location 

d = distance 


R 


m 


Lat. corr. 206265 

6378 Km. (mean radius of earth) 


long. Corr AL" =d_Sing^ 2 06265 
These mapped ground water well points were then located 

on a ground water board map of the region (scale 1 " - 4 mile). Those 
points which feel on any of the six lines (for microdensitometer study) 
or in the near vicinity, were chosen for calibration. The light density 
at these chosen points were ascertained from the micro-densitometer plots, 
for the Gujrat region (lines 4, 5, 6). Only the Gujrat region was 
considered ' as the CCT study had earlier indicated a 1 jack of appreciable 
correlation in the Rajasthan region (see Digs. 2*»5)« The four plots 
were made on ordinary (see fig. No.11-l4) graph paper by plotting ^1^^| 
density versus depth. The first three plots were for b±| 


KLn. 


A 


8(951 
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Table — 8 

Field data for Gujrat region 
(Lines 5, 6 Densitometer method) 


SI* 

No. 


LINE 

NO. 4 

J LINE 

NO. 5 

~T 

LINE 

NO. 6 

X 

J Density 

I 

( Depth 

i Cm.) 

£ . 
j Density, 

5 

f Depth (m. ) 

i 

nr 

i 

5 

Density 

j 

| Depth (m. ) 

1. 


0.796 

5,60 

0.745 

5.90 


0,921 

5,40 

2. 


0,699 

6.50 

0,639 

6.50 


0.903 

5.50- 

3. 


0.620 

0.00 

0.620 

7.20 


0.757 

5.60 

4. 


0.629 

0.60 

0.670 

7.00 


0.796 

5,70 

5. 


0.602 

9.30 

0,602 

0.90 


0.605 

6.00 

6. 


0.620 

10.00 

0,605 

9.20 


0.670 

6.30 

7. 


0.620 

10,00 

0.602 

9.50 


0.630 

7,10 

0, 


0.620 

10.00 

0.620 

0.50 


0.640 

7,60 

9. 


0.629 

9,00 

0,620 

7.60 


0.630 

0,50 

10. 


0.620 

9.40 

0,650 

6.90 




11. 


0.630 

9.00 






12. 


0.670 

7.70 
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individually and the fourth plot was for the combined data of lines 
4> 5, 6. Each of these four graphs established that at water depths 
deeper than 10 m there is so much ’scatter' on each graph that it can 
be concluded that there is no appreciable correlation beyond this depth, 
between the light density (satellite reading) and the water depth* Each 
graph showed an exponentially decaying curve of this mathematical forms 

"h.y -j- r» 

y = e . ' . The equation of each curve was 

established by finding the values of 'b' and ’ c' for each curve by the 
method of least squares curve fitting. Please refer to (Appendix - III ). 

The calibration curves for lines 4, 5, 6 and combined plot 

for lines 4, 5, 6 (master plot) are given in Pigs. 11-14 and. the curve 
equations and correlations coefficient (see appendix 3, are given in 
table 9.- 


Table No. 9 


line-! 


E 


O 

,nel 


Calibration Curve Equation 

B,C 

b~ r — 


0.04761652 
0.04636945 
0. 10596548 


0.02515450 

0.09178096 


{ + 0.36929877 


0.05901143 ! + 0.05047683 


Coefficient of Correlation 


- 0.88415922 

- 0.78772666 

- 0.79501882 

- 0.79258990 


Calibration Curve equation, y = e 


b x +c 


3.1.4 Comments s 

i) The negative values of the coefficients of correlation .> 
indicate that ’y’ decreases with increase in ’x’. 
ii) Possible explanation for the light density to decrease 

with water depth is a.n exponentially decaying curve may be 
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as follows! water does not reflect radiation in Band — 7, 
hence a surface water body such as a lake, in Band — 7 would 
theoritically appear opaque on a positive film ( i,e., on a 
’positive 1 negative). Its transmittivity, T, would be 1 O' . and, 
therefore, its density, D, would be infinite according to 
Eq. (1). As the ground water depth increases, the effect of 
soil moisture decreases and the same soil shows greater 
reflectance and, : therefore, a greater transmittivity on a 
positive film, and hence a decreasing density. Hence a 
micro— densitometer plot of any line on a positive film of a 
particular area, in our case the Gujrat region, would show a 
high density value at very shallow ground water depth, which 
decreases exponentially with increase in the depth. 

iii ) The high value of '-r* indicates an excellent correlation 

between the light density and the water depth upto 10 metres. 

l 

The plots for lines 4, 5 6 can be treated as auxilliary 
curves or verification curves for checking the water 
depths given by the master plot of lines 4, 5 6, though 
the latter by itself has a very high coefficient of 
correlation ( —0.793 ). ( The surface area covered 

by the master plot is approximately 185 Km x 90 Km. of 
Gujrat region). ( See Figs. 11 to 14.) 
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3.2 Prediction. Part (Micro-densitometer approach) t 

Only the Gujrat region phe given imagery was considered for the 
prediction part because the given Rajasthan region indicated no appreciable 
correlation during the calibration phase .(see 'CCT method* page 9 ). 

In the CCT prediction process' the given Gujrat region was 
subdivided into Regions (l) (2) & (3) (see page 19 ). In the densitometer 
method the Region (2) & Region (3) { t or lack of contour points in the 
given depth contour map) were combined as Region (2), & it was demonstrated 
that the satellite survey is not reliable for such regions (see' Results ; 
densitometer method'). 

A number of points were taken on a 1 " = 4 miles scale map of 
the Gujrat region supplied by the local ground water board authorities. 

This map had ground water contours plotted on it by conventional methods. 
Next each of these points were alio ted to one of the three lines (lines 4, 

5, 6) and the position of each point with respect to the corresponding 
line noted. Also the ground water depth at each point was recorded from 
this map, if necessary by interpolation. 

Next the relevant line corresponding to each point on the 
map was considered and the corresponding light density value was noted 
from the microdensitometer plots-. Then this density value was utilized to 
predict the corresponding ground water depth from the master *alibration 
curve as well as from the auxiliary calibration curve for each line. 

The ground water depth(s) so obtained were then checked by comparison 
with the depths obtained from the conventional ground waster contour map 
(see -Rig. 15). This completed the prediction of the ground water depth 

for one point on the map. This prediction process was repeated for the 
other points taken on the map for., ground water prediction^ The mean errtirs 
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Table -10 

Field data and calculation for prediction : 
(Gujrat region 1 ) 


SI. j ~~1 I I I Pred- I Inter -r 

No. $ Longituqe J Latitude {Line {Density { icted 5 n 

{ (Degrees) {(Degrees) 5 No. { 5 depth 8 a ? ~ n \ 

5 { { { { line 456 { Ciru ' 

I I | L |m. ) , P { 


Deviation 

|P 


; x 


(m. ) 


■Di 


REGION 1 


1. 

71.014 

24.543 

4 

0.638 

8.47 

9.00 

0,53 

2. 

71.840 

24.520 

4 

0.593 

9.71 

9.27 

0.44 

3. 

71.064 

24.518 

4 

0.602 

9.46 

9,54 

0.08 

4. 

71.884 

24.514 

4 

0.595 

9.65 

9.82 

0,17 

5. 

71.542 

24.291 

5 

0.678 

7.45 

5.00 

2.45 

6. 

71.579 

24.284 

5 

0*678 

7.45 

7.00 

0.45 

7. 

71.626 

24.277 

5 

0.638 

8.47 

7.00 

1.47 

0. 

71 1 765 

24.253 

5 

0.619 

8.96 

7.00 

1.96 

•9. 

71.777 

24.295 

5 

0.602 

9.45 

9.00 

0.45 

10. 

71.412 

24.044 

6 

0.638 

8.46 

7.00 

1.46 

11. 

71.444 

24.042 

6 

0.619 

8.96 

9.00 

0.04 

12. 

71.535 

24.026 

6 

0*638 

8.46 

9.00 

0.45 

13. 

71.602 

24.014 

6 

0;745 

5.84 

7.00 

1.26 

14. 

71.693 

24.000 

6 

0.745 

5.84 

7.00 

1.26 

15. 

71.793 

23.984 

6 

0.553 

10.90 

5.00 

5.90 

IS. 

71.884 

23.969 

6 

0,638 

8.46 

7.00 

1.46 

17. 

-71.479 

24.040 

6 

0*699 

6.92 

8.00 

1.08 

18. 

71.488 

24.033 

6 

0,678 

7.45 

9*00 

1.55 

19. 

71.576 

24.019 

6 

0.854 

3.53 

8.00 

4,47 

20* 

71.650 

24.007 

6 

0.678 

7.45 

6.00 

1.45 

21 . 

71.743 

24.991 

6 

0.699 

6.92 

5.00 

0.08 

22. 

71.843 

23.976 

6 

0.745 

5.84 

6.20 

0.46 

23. 

71.600 

24.014 

6 

0.769 

5.31 

7.00 

1*69 

24. 

71.000 

23.980 

6 

0.699 

6.92 

5.00 

1.92 



Table — 11 
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Field data calculations for prediction 2 
( Gujrat region 2 ) 


Si. 

NO. 

1 : — : 

» Longitude 
| (Degrees) 

i 

4 

5 Latitude 

} (Degrees) 

1 5 

1 Line 
{No. 

1 Density 
| 

1 Predicted 
J depth 

5 line 456 

J (m. )P 

mSH 

Hilne 

1 — : — : — 

« Deviation 

Jx=|p-d1 
j (ra ? } 

REGION 2 







1. 

72.000 

24.512 

4 

0.678 

7.54 

15.00 

7.46 

2. 

71.992 

24.514 

4 

0.602 

9,46 

13,00 

3.54 

3. 

71.983 

24.517 

4 

0.600 

9.51 

11.00 

1.49 

4. 

72.005 

24.514 

4 

0.602 

9.46 

17.00 

7*54 

5. 

72.011 

24.512 

4 

0.638 

8.47 

19.00 ' 

10.53 

6. 

71.819 

24.249 

5 

0.709 

4,95 

11.00 

6,05 

7. 

71.851 

24.239 

5 

0.658 

7.96 

13.00 

5,04 

0. 



5 

0.678 

7.45 

15,00 

7.55 


71.893 

24.235 






9# 



5 

0.678 

7.45 

17,00 

9.55 


71.954 

24.228 






10. 

71.977 

24.221 

5 

0.638 

8, 47 r 

19,00 

10,53 

m 

72.093 

24.200 

5 

0.638 

8.47 

19.00 

10.53 

12. 

72.219 

24.179 

5 

0,553 

10.89 

17.00 

6.11 

13. 

72.263 

24.172 

5 

0.769 

5.29 

17.00 

11.71 

14. 

72.3Q2 

24.165 

5 

0.658 

7.96 

15.00 

7.04 

15,. 

72.321 

24.184 

5 

0,658 

7.96 

13.00 

5.04 

16. 

72.339 

24.160 

5 

0,619 

8.97 

11,00 

2.03 

17. 

72.358 

24.158 

5 

0.678 

7.45 

9.00 

1.55 

18. 

72.388 

24.149 

5 

0.585 

9.94 

9.00 

0.94 

19. 

72.430 

24.147 

5 

0.585 

9.94 

11.00 

1.06 
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and 'fee standard errors were computed for Regions Cl) and ( 2 ) to assess 
t e accuracy of predict ion in each rngtonC see cage 19)* 

Once it was established that this Method provided a good aeournt 
of predict long the need of 001:3 oaring with each predicted ground water 
depth with the conventional ground water contour sap of an area is no 
longer necessary and the calibration prediction ' curves can he reliably 
used to predict ( or update old convent tonal maps) shallow ground water 
depths in umtaopedfereas* 

aroa ard i3i: imr c ?iic.io*s:-rsi?q.r v r^ a?? ioaoi ) t 

Region(l) 1 ( Sec table - 10 ) 

Huatber of points Sfe«24, 

Standard error , (f m 1.38 m 

Hie 3 - (T test oh owsthat we lave to reject any ac>4*14 m 
hence w© reject the two observations c ocresoondimg to serld. nos 15 
and 19 in taMe **10* Therefore, we have t 

IT * 22, (T» t 0*71 tt 

v*e get a good prediction accuracy, therefore, in Region Cl) a® 
the standard error t»£ 1,00 n* Region Cl) is mostly Alluvial with 
ground water depths less than 10 m, where our satellite survey method a 
are most effective# 

M&JmJM lL» ( ****• ) 

IT « 19, <£» 3*56 a 

The 3 - <T toot shows that we hare to reject any x > IQ *68 a, hence cm®, 
observe ion corre® pending to serial no* 13 in tala® ~ 11 is discarded# 
Therefore, wo got s £ « 18, 0" * 3*39 a# 

The high value of the standard error indicate this; the satellit 

survey ©ethdds are unreliable la Regioa(2) # which is mostly hilly 
and with water depths greater than 10 m# 
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3,3 Sources of Error ( CCT and Densitometer methods ) T 

1, The normal water level in some of the wells in the Gujrat 
region was lowered by Q,5 m, — 1,0 m, due to pumping 
operations, when the field data was collected. Hence the 
water levels in these wells represented a lesser values of 
the actual ground water depths. 

2, Ideally, the satellite survey and the collection of 
field data ( as well as preparing a conventional water 
depth contour map ) should be done simultaneously. In 
practice,, it is preferable to collect both types of data 
as close to each other as possible. Otherwise, the 
effect of seasonal variation of ground water depths 
comes into play, and the information recorded by a 
satellite and by a field observation, at any particular 
sport, will not indicate the same depth, but different 
depths due to seasonal variation. In this study, the 
satellite survey was done in June 1977, the field data 
was collected in December 1979 and February 1900, and 
the conventional water depth contour map was prepared 

in December, 1970, The correction due to seasonal 
variation could not be applied to the water depths due to 
lack of data. The low value of the standard error achieved 
(for ex,,+ 0,99 m, ; CCT method ), without correction 

due to thy seasonal variation, indicates the high 
accuracy possible once the seasonal variation factor is 
taken into account, 

3, To pin point the latitudes and longitudes, a 1" = 4 miles 
ground water contour map was used, for the CCT method, 

A higher accuracy could have been achieved with a larger 
scale map. 



C3*APTK * 4 

Disc^ssiOH* coic:r,usiog a HD B^ortr^cariQi® 

4*1 disc ssioir. ooacL r 3icm aid r^aw^imiom t 

H m well Ota ground water depth be estimated from satellite 

data? 

For -fee given 3a j as than region, the application of the GOT 
met od yielded no appreciable correlation between Ihe ground water 
depth and the satellite data, using the sent tot plsasl method and the 
linear averaging method ( 3, 5 7 niseis array), This indicated tha- 
ws oan not reliably estimate water depths from satellite detain this 
region, which is mostly sandy and with water depths in excess of 10 m * 
The hilly terrain In the lower la last ton region, i,e* Jlount Abu hills, 
were not Included in this study because those prediction methods are 
not reliable in hilly areas. 

Best we studied the Gujrat region given in the imagery with 
the help of the CCT method and the micro-densitometer method* fhree 
different types of C for legion© 03, (2) f (3), see page 19 3 wore 
classified depending upon terrain C alluvial or hilly 3 and mater 
depth { greater than or less than 10 a )» The oallbraion prediction 
yielded very promising prediction accuracy { standard ©rr©r*0*99 ml, 
in case of OCT method andstandard error 0*71 a in case of densitometer 
method) C sc© page® 20, 40 ) in Series (l), where the terrain was 
mostly alluvial with water depths less than 10 a* It wasin this type 
of region that we had the best results* On the other hand, as antisl- 
patsd, both met" ods yielded poor results ( see sages 21 y 40 ) for 
Hs/;ion* (2) and ( 3 ), where the terrain was mostly hilly and water 
depths greater than 10 a* 
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A study of field data { see p* 29 ) is the Gujrat region, 
of dep-fhs greater than 10 m, established that the CCT method yields 
a poor accuracy (standard error » 2*63 m, see p«20 )• 

Thus we can conclude that both the CCT method and the 
micro -densitometer met od are best applied to alluvial soil with 
shallow water depths not exceeding 10 m. A few points of interest 
are discussed here* 

CCT Method Vis-a-vis the ^icro-densitometer Method * 

For the given Gu e lrat region* the correlation coefficients 
and the standard errors - 0*939 m and 0*99 m ( legion (1) ) iscase 
of the COT met od, and -0*793 m and 0*71 n* in case of micro-densito- 
meter approach* It was found that ground water depths deeper than 
10*0 a cannot be estimated because of a lack of correlation with the 
satellite data* Two points are to be noted* First isthai the above 
values of the correlation eoeffici ents { 0*8 - 0.9 ) indicate a 
s strong correlation and the second point is that it® value is higher 
for the COT approach* (The correlation coefficients varies f rom-1 ) to 
+1 and the closer to unity it is, the stronger the correlation.) 

A correlation eoeffici ent of 0-0*5 ( test at $$ level of confidence) 
indicates no correlation to a negligibly small correlation while 
0*7 - 1*0 indicate a good to an excellent correlation). Again the 
standard errors are very reasonable considering the soirees of 
errors ( see ’Sources of error* )* The fact that depths could he 
predicted to a metre (CCT approach), pointsto the higher degree of 
precision possible with the minimization and removal of these sources 
of errors .The standard error is lesser for the CCT approach* 

The reason for this is that in CCT method, we have for more d*ta 
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available ( collected in lines and pixels ) about every 80 m we have 
a reflectance reading - while in the microdensitometer method such 
precision was not available in the calibration part stage. The GOT 
approach is superior because it offers an enormously large wealth 
of data. ( 1 line about a 185 Sh ) contains 2656 pixels, each pixel 
about 80 m square ). The next LANDSAT has a resolution of about 
30 m and a ^French land use satellite is being planned for a resolution 
of 10 to 20 m.with th e s e. With these satellites, the CGT method will 
offer an even more accurate picture of the ground water profile and 
an accuracy of within half a metre in area^bonducf ive to satellite 
survey can be realized. This is not to belittle the microdensitometer 
method, which by itself is good enough to give good estimates but 
only that a higher precision can be obtained with the C®T method. 

Central pixel method vis-a-vis the Grid averaging methods (OCT approach): 


In the CCT method, in the study of Gujrat region, the 
correlation coefficient decreases from the central pixel method to 
0.55 for the 7 pixels Grid pixel method and hence the degree of accuracy 
gets lesser. The reason for this trend is that we are assuming the 

ground water depth ( which remains constant over a limited area only) 

to remains at a constant level over an increasing area as we move 

from the central pixel method to the 3 pixels, 5 pixels and 7 pixels 

Grid averaging. The precision therefore decreases. In the 7 pixels 

land case, the correlation coefficient is only 0.55, indicating a 

very poor correlation. Thismeans that the ground water depths should 

not be averaged beyond a ground surface distance of about 400 m allowed 
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the 5 pixels grid method (each pixel 80 m, length of a 5 pixels grid 
5x80 = 400m). Otherwise there will result a negligibly small 
correlation or no correlation as in ca.se of a 7 pixel grid averaging of the 
Gujrat area. Also it means that the central pixel method is most 
accurate in the COT approach.. The restriction of 400 m. on the grid 
averaging also means that a remote sensing satellite should have atleast 
a pixel resolution of 400 m to map shallow ground water depths. The 
first Indian Remote Sensing Satellite, Bhaskara - 1, has a. pixel 
resolution much larger hence lacks the precision required'. It is 
hoped that with the Indian Space Research program maintaining its 
ambitions program, I8R0 will in a not too distant future put a 
Bhaskara series satellite of pixel resolution well below 100 m. When 
this happens..^ indigene ous satellite survey based on the OCT and 
microdensitometer methods explained in this study. will be realized. 

Till we achieve the required pixel resolution in our very own Bhaskara. 
etc., ISRO can procure LANDSAT data to develop a satellite ground 
weter depth contour map of most parts of India. 

Satellite Survey vis-a-vis Airc r aft Sur vey ; 

Aircraft Survey ^Versus Satellite Survey 

The objective of ground water depth prediction 

can be achieved using both an aircraft and the satellite with MSS 
equipment. But the advantage of a satellite is the larger area coverage 
(185 Km x 185 Km”. LANDSAT) compared to an aircraft, A satellite survey ; 
more suitable in a ground water depth estimation in large areas, 
which would require a large number of flights (and camera exposures) 
by aircraft. Aircraft survey is suited for covering small local 
areas or for those projects where steroscopic photographs are needed.- 



It- is obvious that to prepare ground water maps of large areas* 
the superior option is a salellite borne survey. 

Summing it up, this study has proved that it is 
possible to predict shallow ground water depths from remote sending 
satellites with a good accuracy. Ground water depths upto 10 m. depth 
can be estimated while rocky areas or built up areas or areas 
under dense forest cover etc, cannot be covered by a satellite 
surveys. The given G-ujrat region was probed with two techniques, 
one using CCTs and the other a micro-densitometer. Upto 10 m. 
depths could be estimated with a high correlation and very reasonable 
standarpd errors especially when quite a few sources of error existed 
Undoubtably the standard error of + 0,99 approach) can be 

refined considerably with the removal of these sources of erro.r. 

The highly encouraging results obtained in' this project highlight 
the tremendous scope of predicting s ^ a -^ ow ground water depths by 
satellite surveys . 



'AWcjOIX - 1 


1 


U>00 

90S 


5 ;,fS pRUGHAh CONVERTS THE LATITUDE AND LONGITUDE OF A PljT NT flrt fHf 
T-CtCO IT! ! is CuF^ESPONuIriO M*t£-W 0 i* 8 l(.R AwD PIXEL DUMBER if' Tbl 
' Ai.EkY A |* T E R CORRECTING FOR THE CURVATURE UF TriE EARTH, 

TOPbiCll DOUBLE PRECISION (A~H*0~Z) , ^ , „ , „ , 

dimension plqoo) ,pphi(200) »xp(20o ) # yp(2oo) , u(2oo) , v( 200 ) 

9, \f ii C 200 > 

OPEN ( UNIT»7 , DEVICES € DSK ' ) 

REAU(7,*)N 

Rc » U i 7 , 1 ) ( PI. 1 1 ) , PPHI ( 1 } , |*1 ,JN ) 

FORMAT (1 4) 

F0RHAT(2FI1 .6) 

DO 905 1*1, L 

0PHI(I)aPPKlCI)N3.14159i/180. 

C A h 1, V A L C KPH T , R , OpHI , OL , P M , RN ) 

D{ t *FL(T)»Ur. 

DpHI = yPHI(I)-OPiiI 
Csnb*(Si!»cnnui) 

S.y, = RM* (DPHt ( 1 WiPHT ) 

U « 1 I 5$Wt* 3 

dai 2=6HR**/) : 

0 ;t=DF1 1/C 'J? 
n ; l 1 s G ; . 1 1 4 . 

•' ?i?=Si *'(OFM )/COS(OHFn 
■■o . 22 = 24* (hM 
O' ?*Dw2i 1 *». 7 1 7./VH72 
f)/sSM tD'-sl *»DY? 

R H 1 'risk Pi 1 1 -0. 

C = C * 3 . t 4 1 5 9 3 / t ft 0 
, SlCaSlN(C) 

Xp( i )sRPDM*MC 

Yl'C I )=RRH I - ( RPR I-DM ) *C n S f C ) 

TiPr.;»,XPCi),vp(xi 

ft I --0D -0021 4750 
rtl=-o ,0i 33 1 120 
C 1=3839.33110000 
D t=0 ,01 43782ft 
m gl*«&*O03f2949 ;.t-_ 

Fl**72 . 1 9580300 / 

0(1 )«A1 

V(I) = Dl+XP(ntFl»IPCI.l+Fi , % 

PRTNT1000,T,PUa) ,PPHI(I),U(n,VU) 

fnuu » 5* f iv r . iy ' f nKIflT wififfa " • . ff t I 


m^Ulkkl’ ■ .KSIP’PKtt'SMdl*: KUTS" 

COOT l RUE 
STOP 

SU8R0 Of I H E VALtRPHI ,F , OPHX , OL , RM , RN ) 

■T«6. J7730J?4 

F*1 /foQ , 8 

DPMI *23. *3. 141593/ 180. 

B#E**2 

SaSINCOPHI) 

Sii*S**<s 

SslISs ” 

Sir* 1 

R* “ 

RP 

rctuIn 
END' C T 

> p R r £ -W'. 1 


f , Ft 1.6, 
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PROGRAM TO READ THE REFLECTANCE VaLi,e.> [H jx^ •!„*>,• r-r 

r ',■■■'■- « - ( ">..*?£! ’ THEM TO THEIR DECIMAL. VALUES » 

integer . > ,ht 2 ,bt 3 ,bt 4 .bt 5 

..- 1 • ’ •!. ^,LA5IC30) ,VL(30) 

^PEN ( UN 1 Tx?o , DEVTCEs ' MTAO * .MODES 'DUMP' 

'; ‘ STZE=R40, DENSITY* 800*) 

OPEN ( U»IT®2 1 , FILES* pr.R? l, DAT' ) 

DATA I3/5H / 

)AT : <ST/30*999/ 

ACCEPT*# «MH,M 
FORMAT | 5 X # ' t* 

PRINT 


500 

1 12 
RR 


t 1 t 

4^4 

R 8 R 


■D.i 
ft b 6 


IS* , 15} 


THE LINE NUMBER OF THIS RECORD 

.1 ,MMW 
NlsM+ 1 ? 

XATsO 

MR EC sO 

^CiP 5 Ns 1 , fi »j f) 

T A ( N ) * I 8 
READ(20)TA 

MQ*-3 

*-Vl 1 

IF ( K ,GT „ 1 } GO TO «u 

TsO 

Tsr + i 

KC~0 

rP’(lA(I)*£0.TB)(3O T U tip 
I f ( I A ( I ) , GE . P ) GO TO 25 

tA(l)sTAfl). ARP. *37 7777777777 

KC® 1 

8n = TACl)/24l28 

BC?*CIACT)*DTJ *2**281/2**20 

8T3*CIAU)-t3Ti*2**2ri-HT2*2**2u )/?**12 

3T4»(IA ( I) -Bt l *2**?8 «hT2*2**?0*8?3*2**12) /2**4 

B 15*11 At 7 )-B" r l*2**?B-BT2*2**2Q-ST3*2**t 2 -814*2**4) 

T F ( 8T5 . GT, 7 ) RT5a- U3T5-R } 

IFC KC #EQ « 1 ) BT'j. »8T1 *1 28 
NJsNO+4 

IFCM.GE.MI) GO TO 222 , 

IF(NO,NE.M)GO TO 222 
XRt TE ? 21 # 500 ) NO # R§1 # RTt »BT3 # BT4 » 8 T5 
FORM AT ( 6 15 ) 

M*»M*4 
GO TO 222 


NRTTEC21,500)UAST 
CLOSECUHJT 


y&§§ 

S'r': 


w * * * d t ) : ■ ■ 

fKAfPK AT+5 

RgAD('2l,500,EflDsf888) (VL(I8) #T8*KAT# JKAT) 
If A' ; 9 T ’• ■ 

TFCKAT.LT* 30) GO TO 77? 

XAT*0 

PRINTAAi # VL 
FORM At (3X (i 30 14 3 

lliii ill Me 

50 I ) 99 

C LOSE ( BN T 1*20 , DEVICE* *HTAO ' #MUOL‘* 'BUMP ' 
9, RECORD SIfcE*84Q. DENSITY* '800') 

FORMAT C i NO OF RECORDS READ* ',15) 

m* 


Vi 'f HI 
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' Prt’VJRAM TO FIND CURVE EQUATION 

Di MENS I UN X(4u},YC8O),A(80) 

NilNseO 

fisTS ■ ' ' 

RBAuUl ,*J fxcn , Tst ,W),(Y(U),Jsi ,N) 

PfifWT 44. N 

H FORHATC15X, 'INPUT OAtA FOR C AL 1 8 R AT 1 0 N " PR ED I C T T ON CURVES*. 

?5X , *S* f 12 ) 

DU 4? T= 1 ,75 

roRrtATUQX, *UI) = ? ,Ft5.8,5X, * Y C 10 = * , FI 5 . 8 ) 

42 CONTINUE 

fi ss N • 4 4 

** CONTINUE 

C»D*0 

ANK*i*I«0+l ■■'■'■ 

03 \ i * M , N 
7,sY(T) 

ACT ) = AT,OG(7.) 

c-Etun 

i CCHTTNUE 




0 = 0,0 
TA» 0 .U 

Ran 5 

00 2 J*M,« 

’4*X( J) 

D0b$* w 
0*n +np 
8»B+N 
AA«W*W 
TA*TA+AA 

Continue 

irOfUtftGE.OGuTO II 


irtMW s LB,3U.HuM*4S7 

P.lT.pr 40 0, iij I •:■■.■■:■■ 

FUPrfA'n (////, tuX f *MTcr(j“UENSITOMETER METHOD* , 1 OX, 

HiS 1 ' NUMBER ( 5 ) * * » 1 2 ) 

j t ? ,i ■■ . ■ ; v ; . 

Cl ■ rt iUE " : ■ 4 ■ 

| iiriMf2 

’ -■•.U , 4 !. ; 

FORMATf//// t luX, f MAGNETIC TAPE METHOD* , 1 OX, 'NUMBER OF 

T .’'T : - 4 ,, ' :444, . ' 4 

BBSCB*C-A#IF#D)/(B**2“TA*ANP) 


rnlfl t tUVf nOfVvliivH 

100 FDR 'tAT ( / // , OX , S WB= * , FJ.*> . R , 7X, 'CC* , FI 5 » 8 , 2X, * NUN* * » 12 } 


ilill 

Yl*EXPff| 

YD J F*Y ( 1 5 **YY 


1 i 0 


.■ V . ,‘W>' •: 


CONTINUE 
C S*O,0 

no do n 

CONTINUE 1 J 


*XDIF«* , Ft 0,4, 

A) 


SlAt ' ■' 


' : 1 
I II '44, 

44 ,4 4 ‘ . 4.4 •: ;4:4,-i>' 44,4- ':\4 ' 4 1 4; 4'. , 


, 1 > , 4 v . 4 # 4 .. 4 !’ 0 : 444 : 4 ‘ : ->S • \ 



• r . " V . . t , ■ ,, 

FJR^TC//,5a, VftVX*',Fl5.R,2X, 'ft VY=',F15.«,2X, 'NUN* 4 

7, S ' ! j4 ss U * a 

XSUH*o*0 ' 

no 50 IsM.w 

0»X ( T ) •'ftVX ' ■ 

R s 0 * * 2 
S*M T)~AV¥ 

Va &**2 

0 ^ 0 , X ( I > » C 1 1 

FUP«AX(/,5X,- f X(U* # ,Fl0.4,2X f 'Y (I )s ' , F6 • 3 , 2X , 'RX=' , 
42*,'gY*Sn?.8,2X,'A.gXY=',F{i> J) 

XSU*J*XSUM+R 

Y3'M«YSUM+Y 

ZSUfS*2SUM+Ag 

CONTINUE 

TF(J,LT,30)J*J-1 

AX*XSUMA1 

A f*YSUH/J . 

XXSTDsSORTUY) 

Y (5TO*s5PT(AYj 
TT3XSUM*YSUM 
S5»SQHT( TT) 

C:jRK*ZSUM/S5 

P6T.« r 600 , XAST n , YYS in , roKR , J , XSuM , YSIlMiZSIJM 
F jRM AT ( / / , 5X , ' X t T ) r.xns ' , po . J , 2 X , ' H I) STL) = V , F6 . 3 , 2X , 
7 'CORRELATION CgEFFIC I Frits' .F15 .8 . / , 1 OX , ' J» ' , 12 , 2X ( 
8F15.8,2X, 'XSUMs',Fl5.R,2X, * XSUM* 1 , fl& . H ) 

♦HINPiUN+1 

GO TO (6,7,8,11 ,12,13, M, 15), NUn 

CONTINUE 

N 3S 1. 2 

GO TO 5 
CONTINUE 
Ms 13 
?U#22 

.continue 

M*23 
N»3t 
GO TO *5 
CafflNUE 
M «|2 
N*42 
GO TO 5 
CONTINUE 
Ms 43 
N*53 
GO TO 5 
CONTINUE 
Mss’? 4 

M#6 4 , 

GO TO 5 
CONTINUE' 

M*65 
Ns?S 
GO TC 


, T2 ) 


FI 5.8 


' XSUM 




x.c|) = 
nil* 
nu s 

mi* 

nus 

X C X ) a 

x(D 8 


P'.eoOOoooo 

5 y"o^o 6 o 6 co 

l#: 8 S 88888 g 

11:88888888 

ililliiilio 


9;90000000 

4.80000000 
(>1800000*^0 
5,90000000 

8 .80000000 
7.^0000000 
5,40000000 

l 8 :?SS 8 SS 8 S 

|.80ig0||| 

Itfoollloo 

iMMS 

R *70000000 
0.90000000 
4.80000000 

fMl 38 

5,00000000 

7,|^||0g.g 


0000 


Aims 9 l) t x " 4 • 


NH.IT o kSA N* 

5.60000000 
6.50000000 

5.00000000 

8.60000000 
9^30000000 

10,00060000 
X 0.00000000 
10.00000000 

9 .80000000 
9 . 4OQy0000 

9.00000000 
7.70000000 

5.90000000 

6.50000000 

7.20000000 

7,800 >00 )0 

8.90000000 

9.20000000 

9 .50000000 

8.50 0 00000 

7 , 60000000 

6 . 90000000 
5, 10900000 
5, 50009001 
$*60000000 
«. 70OOOOOO 

6.00000000 
5,30909000 
7.10000000 
7.80009000 
8 .80900000 
9 , OuOOOOOO 

4.90000000 

6 .80000000 
5 ,00900000 

3 : 9 smm 

" ,40000090 

10,00000000 


ti ^ 
Yf I)* 

r(I)* 

YU) = 
YU)- 

Y ( JL ) * 

mi= 

YU)* 

m)= 
¥(n = 
urn* 

XU)s 
y a i * 
y c i ) » 
y (il* 

Y ( X ) = 

StliS 

UV‘ 


9.79600000 

3:13333333 

3:33133333 
0,62000000 
0,62000090 
0.62000009 
0. 62900000 
0. 62000000 
„• 000 
0.678000U0 
0.74500000 
0.63900000 
0.62000000 
0,67800000 
0.60200000 
0.585000u0 
0.60200000 
0 , 62-0090-00 
0.62000000 
o , 6 5 b 0 0 9 0 > 
0,92100000 
0.90300000 
O,757OOOo0 
0,79600000 
0,68600000 
0.67800000 
0.63800000 
U. 64800000 
0,63800000 

58.00000000 

19.00000000 

27,00f 

19.001 
‘ ' 3€M 



M ! C JU - tU c> f TO fiETER METHOD LIME NUH&ER C S )*457 


r} H s — o . U S 9 

01143 CC* 0.0*>a47683 

. WUN* 0 


x c n = 

5 » 6 1) 0 0 

XX* " 4,7217 

X L> T F * • 

0,8783 

U I J * 

. , nu 

Yi * y .7558 

YD IF* 

0.0402 

X C I ) * 

6.50 0 0 

XX* . 6.9239 

XDIF* 

*0,4238 


' 0,6590 

YY* 0,7167 

YDtF* 

-0.0177 

■ C I ) * 

■ . u-- : : 

XX* 8,9561 

XOIF* 

*0,y5bt 

Y l T } = 

0 . 6 2 0 0 

Y Y* O.oStjO 

YDTF* 

*0.0360 

X ( T j = 

8 ,60o0 

XX* 8,7119 

XDIF* 

-0.1119 

v ( t ) - 

0 . t> 2 9 0 

YY* 0.6332 

Y 0 1 F * 

-0,0042 

xtn = 

9, JO 00 

XX* 9, 455 ^ 

XDIF* 

-0.155 4 

Y CD* 

0,6020 

YY* O,t>07 5 

YD IF* 

-0,0055 

x a j * 

i-j.yooo 

XX* 8.9661 

XDIF* 

1,0439 

YU)* 

0.6260 

YY* 0 .58 JO 

YDTF* 

0,0370 

X ( T ) * 

10.0000 

XX* 6 ,95b 1 

XDIF* 

1,0439 

YU)* 

0,6200 

YY* 0,5830 

YD IF* 

. . ; 37 0 

XU)* 

10.0000 

XX* 6,9561 

XDIF* 

1.0439 

YU)* 

.»200 

YY* 0,5830 

YDIF* 

O.O370 

m 
■■"*** ■ 

**# : 

If 

9 . W 0 0 0 

XX* 8,7119 

XDTK* 

1,0881 

. YU)* 

*.0 290, 

YY* 0,5899 

YDIF* 

o, ;i 3y i 

’ «!)* 

9, >0 

XX* 8.95oJ 

X D 1 8 = 

M. 5 YU 

■ -u 

' YU)* 

0,6200 

YY* 0,6040 

v. ; tv = 

0.0160 

XU)* 

9,0000 

XX* 8 .474 1 

XDIF* 

0.5289 

YU)* 

0,6380 

YY* .0,6184 

YDIF* 

0.0196 

XU)* 

7.7000 

XX* 7.446? 

XDIF* 

0.2593 

VCD* 

. ? 

YY* • < .<’577. 

YDIF* 

* f- '36*. Ah'.0 2 jw 

0.01O3 

X u J * 

b.yooo 

XX* 5. *437 

* YY* 0.7425 

tD IF* 

0.0563 

YU)* 

0,7480 

YDIF* 

wto'f.4 - ;* • 

0.0025 

XU)* 

6,5000' 

XX* 8.4446 X0ir* 

U 2,..,. .7>'. ' ..’..> -4 >..r' - . ...... 

-1 .94*36 

YU)* 

.'.0« 70 

. ’ ' ii 

YY* 0*716? 

YDT.F* 

-0.0777 

X £ T ) * 

/.2000 

XX* 8*9561 

XDIF* 

-1*7561 

YU).* 

. v.2 .* '■ 

Y Y * , . - - 7 7 

YDIF* 

-O.0677 

XU)* 

7 . s U-> 

XX* ' ' 7,940? 

XDTF* 

0.3593 

•YU>* 

- 

■ "YY* 0.6638 

YDTF* , 

0*0142 

X<I>* 

«u##o 

XX* 9*45.54 ' 

XDIF* 

' *0.5554 

YU)* 

o.-ote 

YY* 0*6221 

YDTF* 

-0.0201 


y.„2<H)0 

XX. *,«*»* 

F* 





-1,2666 

-0.0532 

-1.J407 

-0.0472 

-1.3711 

-0.0538 

4VV ,;U 

-0.6076 

-0.0237 
• 0.0289 
0,0011 


K(I)« 

n*)= 


'■ - v¥(!3* 
X(t)* 
-■ ¥(!)» 


8.5000 XX* ■ 8*4711 
0.6380 n* a.0369 


YDIF* 


YSTlf 


113 ,, 

:>** • 4? 


STAr; 


CQRR' 


M I CK0-n« f 18 1 TO Y»ri'CT METHOD 

it .'!*■ : ‘ ■ 'tv-, t-V i i '■ 4 / : S'. ■t.VS VV ' "■! i't ' ste: 


-0.0I5I0MO 

,4.4' 0 3 3 : 

0.?4u6 " 

,.w,- 

. 0.7156 

ia s tsiS* 

9,5110 | 


„ ,-0.047 63' §52 ' I 
'' S.6000' 
yen* .,,0, , ..790 , 0 


Y( T)* 

X ( T ) = 

Y ( X ) * 
XU )* 
Y(X)a 
XU)* : 
Ytt)* 
XCX)* 
a ci)* 

X(I) = 
’((I)* 
XU) = 

Y l T ) * 
XU)* 

v( n* 
x t n * 

Y(I) = 
X(I) = 

ten*' 

X(T>* 


0. 5.550 

- 


•Yi)T,-s 

9.5000 

: = 

lillllfit:-) 

xoiir* 

0.6020 

yy* 


YDIF* 

8.5000 

* V - 

8.9561 

XDIF* 

0,6200 

n* 

0.6369 , 

YDIF* 

7,6000 

XX* 

'0 ,«54i' : 

XDIF* 

0,6200 


■Sl||ff7tt 

YD IF* 

6.9000 

XX* 

7.948) 

XDIF* 

0.6580 

¥¥* 

0,7000 

YOTF* 

6,4000 

XX* 

2.2499 

XOIF* 

0.0210 

r<= 

0.7648 

■ YOIF* 

6.5 0 0 0 

XX* 

2.5 34:4 

XOTF* 

0,9030 

YY= 

0 . 1 6 0 3 

YD IF* 

5,6-000 

XX* 

5.5730 

XDIF* 

0,7570 

. 

0.7558 

YDIF* 

5.7000 

1 XX* 

4.7217 

XDIF* 

0.7960 

fy* 

0.7513 

YDIF* 

6.0000 

XX* 

7,2666 

XDIF* 

0.6850 

V * 

0,7382 

YDIF* 

6.3000 

XX* 

7 .4407 

XOIP* 

W.6780 

• VS 

0.7252 

Y01F* 

7 . 1-6 ofe 

XX* 

8.4711 

" ..'ioir* 

t-'V;-': : : s:vt . 

v .6380 

YY = 

0,6,916"', 

YDIF* 

7,6000 

XX* 


XOIF* 

y,*480 

YY* 

. ''0.6,717 

, 

YOTf* 


- 0 ,09 F t 

0.0446 

16 

-0 .4501 
-0.01*0 
-1 . 3561 
01 7 

-i.: 6 1 

-0.0420 
3.1501, 
0 . 166 ? 
2.91 56 
0.1427 
0.0270 




Aw. DCV, * 1 ,45 1 YSTA*. DEV.« Q 
TON DOSFf tCIEfoT* -0,68415922 


LINE NUf5DE«<S)« 


MICRO-DENSITOMETER METHOD 


O, 091 780 9$ N.UM* 2 


(1,1*200, ¥Y* 0,6534 YDIf* 

7.8000 XX* 6>4013 X0IF* 

*** 3.-y*s-' v -' T f' 


Y 1 1 ) * 

xu>* 

Y(I) = 
X(I)* 
YU)* 




. . 


v {, T } a 

*0.6200 

11* ■ 

0.6 663 

YD IF* 

-0.0463 

:<(!) = 

8.6000 

XX® 


: >IF= 

-0.6083 

II 

- . 0 

Y¥* 

SSli®: 

Y0XF* 

■yy . .;Y- - : 

-0.0155 

1 


y . j o . > 

XX s 

ID, 4297 

y'lplF*, 

-0.82-7 

Y( T } * 

20 

yy* 

0.6263 

YQIfs 

-0.0243 

xu)* 

10.0000 

■|g;y y -jMpy • • v''r%' y " ' ■ 

XX® 

9.5UO 

XPIF* 


;YU>* 

0,6200 

ys 

0.6<>57 

YOTF* 

■).0 «. 4 1 

X C I ) * 

to .0000 

XX® 

9,5110 

XDIF* 

0.4890 

YU)* 

0.6200 

n* 

0.60.57' 

YD IF* 

0.0143 

X ( T ) * 

i 000 

XX* 

9.5110 

X o T F * 

0,4390 

X 1. 1 ) = 

0,6200 

XT* 

' 0,6057 

YD IF* 

0.0143 

x a j s 

' 9 . H O 0 O 

XX* 

9,2 Ob 1 

> D T r * 

a.5‘U7 

Y(I)s 

v - . 6 2 y 0 

YY* 

.0,61 15 

y oif* 

0.0175 

xu>* 

9.4000 

XX* 

9,5110 

XDIF* 

-o. mo 

YU)* 

■) . .>200 

yy* 

U.6231 

YDIF* 

-0.0033 

XU)* 

■ , '0 00 

XX* 

8.9100 

XDTF* 

0,0900 

YIIJ* 

. , - 180 

YY* 

0,6353 

YDIf* 

O.0O27 

xcn* 

7,7000 

V XX* 

■7,6329 

XDJF* 

0 .0,67 1 

f*t 
#*% ~ 

■w 

If 

0.64:80 

'Y* 

. ■ 0,675$ 

YDIF* 

■ 0.0022 


rll • ,iM ; l> ’ .<■* 

XU)» 

:■ * " ■ ■■ 

XX* 

4,369*3 

XDTF* 

1.&310 

ilui* 

0.7457 

6,5000 

i Iks 

¥¥* 

XX* 

v ■ V' 1M : .■ 

YY* 

0,6939 

7,6790 

YDIF* 

0.051 t 

|J § 1 1 §§ 

YU)* 

0,6390 

o,.>7 99 

Yt.ir= 

' 

-0.,O359 

X ( I ) * 

> . ’ . ' 

XX* 

iilill 

xaiF* 

-1 • 1299 


Bfl 


XU)» 
v II j~ 
X 1 1 .) * 
V(T)* 
XU)* 
>".» 
I*- 


8SIBSI 




9 , 5 9 0 0 
0 • 6 0 2 0 
8.5000 
o ,*.>200 
/ ,40O 0 
y*42 ao 
0,9000 
0,6580 


XX* 

YY* 

XX* 

m 

XXf 

ft* 

xx* 

¥Y* 


8 ,96: 5 § 

i§Ift§ 

; 9 *3399 
0.6151 
8,3299 

0. 6413 

1, "' ■ n 

10:, 6625 


XI) Tp a 
YDlf* 
XDIF* 
YDTF* 
XDIF* 
YDIF* 
X:)If= 
YDIF* 


0 , 534 ? 

4? 

O', 1 ?o * 
0.0040 
- 0.7299 
•0.021 3 
-0,1171. 
-0,0045 


xsn«. *>:?.* 1,205 YSU>. DEV.* 0.071 

ruPWFl.4yT.-i J COEFFICIENT* ->,78772666 


MIC HO- i> £ oITO Pi FTEB METHOD 


Ol«E VUMBEKlS)* & 


-0,10596548 C 

C<U; : " I" '' 

0,36929877 

SUN* 

3 

X n : ) * 

5.4000 

XX* 

4.2617 

XDIF * 

1 . n i 

YU)* 

0,9210 

YY* 

0.8163 

YDIF* 

0 , 1 0 4 7 

xcn* 

5,5000 

y XX* ■ 

4,4480 

XDIF* 

l. H,u 

YU)® 

0,9030 

YY* 

0,8077 

YDIF* 

0-. 0*5 3 

xcp* 

5,6000 

XX® 

6.U23 

XOff* 

- >.6123 

YU)* 

0 . 7 5 7 < 

YY* 

0.7992 

f«F* 

-0.0422 

XU)® 

5,7000 

XX* 

’C: : ::;5l634ffU 

XDIF* 

0.0618 

YU)* 

0,7900 

YY* 

0,7908 

YDIF* 

0,0052 

XU ) * 


XX* 

7,455^ 

XDIF* 

-1.0555 

■■ I.'ux-iu' ;■ ' l! 

YU)* 

0.6850 

YY* 

0, 

YDIF* 

-0,0811 

XU)* 

6*3000 

XX* 

7,1524’' s 

XDtF* 

-0.8524 

Yen* 

0.6780 

YY* 


YDIF* 

-0,0641 

XU)* 

7.1 000 

XX* 

7*7 262 

XDlFrf 

-0,6262 

Y,no* 

0,6384 

H* 

Iliil 

YDIF* 

-0.04J8 

|| ./ 1 , 

$;t u=. 

7.6000 

■ XX* 

7,5795 

XDIF* 

Q* 0205 

i! 

*■* 

3*- 

■ 7,6480 

YY* 

1111 I! 

YDIF* 

0.0 H 1 

ixu)* 

d .5000 

XX* 

7*776 2 

i.X: ! i I 

XOIF* 

b. 4.77 3 8 

tit)* 

• 0,6)80 

' YY* 

ll'ilil 
.. . , 

YOIF* 

O.05O2 

■ 


lOOy.i X'""' ; 


'■ vv 

# 





NUM&EB OF PIYEECS1* ? 


.MAGNETIC TAPE METFOLj 


0, 16010599 C.c» a., 4129-453 
mu* ' '9 t 9O00 , , XX* HI. 39«P 

VCt)* 59,0000 YY* *4.49«f 

fCX)» -1,8000 XX* 4,2353 

Y(I)s 22,0000 xin 24*081? 

XII)* 8.8000 XX® 5.2787 

Y;C?)* ■ 28*0000,, yy* 33*470 

xtn* : a. oft op ' xx* s,5t4%' 

*«)» 5,*M» XX* 


XOlf* -0,496ft 
y0!F* 4*5), 12 
XDlF* 0.5647 
YD If* -2*081? 
XU IF* 1.521 3 


NWWiFf* OF pIXEtCrO* 1 


'AijNFTTC i ft PE N 6 f BOD 


08* 


v . 2 0 

175720 C 

■ s 

-1,92411720 

MUfrl* 4' 


XU)* 

9 *fp,0 

>: > » 

lillllf 

’ ; H's 

-0 .60 8 6 

un* 

56,0000 

, yy* 

50*4768 

' r*-- 

7.5232 

; ,xm« 

* *80,00 



X 0 i i; -® 

-0.26 72' 

' ytn* 

3 9.UOO0 

„ yt* 

'18.0393 

lilll 

. . - ;* , 7 

" yen* 

0,8000 

; S 

to. 7989 

XP1F* 

0.0011 

y ( r ) * 1 

77.0000 

yy* 

27.0062 

YD IF* 

-0.0062 

X C 7 ) * 

9.OO0C 


5.0572 

XDIF® 

-y . v5?2 

ye u* 

i 9,y0x/0 

',yy*' 

16.7821 

YD IF* 

0.2)79 

X(T) = 

5.8000 ’ 

XX* 

5.5532 : 

XOTF® 

0.2468 

V(I)» 

7 :i . c '■ ■ c 

yy= 

22.0720 

YD IF* 

-1.0720 

XU)* 

7. ‘ . r- 

XX* 

6,2351 

XDIF® 

1.4849 

y (i)*- 

74,0000 

yy* 

32.3835 

YD IF® 

— , .v:.*S 

X(I)s 

5. »V<. 1 

XX® 

6,2151 

XDlF* 

-0.8151 

Y t T j * 

74,0000 

¥Y = 

20,3607 

YDIF* 

3.6393 

XU) * 

10.0000 

' XX® 

10.6733 

X D I F* 

-0,6733 

Y ( I ) * 

59.0000 

yy= 

51,5055 

yor» • 

7,4945 

Ul)s 

8.7000 

XX* 

7.9414 

XDIF* 

0.7586 

y<n* 

34*0000 

, yy = 

39 ,6228 

YDIF* 

-5,6228 

CORRe£a)'Jo£ CUFF* l 

cIfJt* 

YSTAN. DEV 
0*9369 

5»g- 4 * 0 




■ : ’ 'gg 

ggggg g, : •;> -.g,, ;;,g g, 



X* 

2.41294530 pi® 5 

, ,xx* 

XDIF* 

YY* 

*4.4S*P YMF* 

g gw : . " 

XX* 

.2353 XOlf* 



v U ) * 

*9,0000 

ftP 

X a ) » 

s , 4 o 0 o 

|| i 

* x c i } * 

3/.0Q00 

yy= 

XU)* 

1 0,0000 

xx* 

yu)» 

57,0000 

y < = 

,x(t)* ;; 

8.rw 

XX* 

YUJ* 

||f| 

43 ,0000 

i¥* 

XSU 1 

cokmclah 

. DEV.* 2.047 

N TOSFFICIENT* 


DAGrlFTT 

C .TAPH- ;-1£ 

fHOD 

0.07443*20 C 

C* 

X ( T) * 


XX* 

Y(!)s 
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APPENDIX-V 


KIE1.0 DAI’ASKAJASXHAN REGION (. CCT METHOD) 

*A 9% .. * ft U -ji* V m I t #% r &«*<•’’ «li » rt m i . # « « 7 <«* *■**, "HT 


S,NU. 

LONGITUDE 

1 DEGREES ) 

latitude 

(DEGREES) 

LINE NO, 

PIXEL 

1 . 

71,610000 

25.350000 

1 / 9 

242 

2 , 

72,25277b 

25,244444 

235 

928 

3, 

72.51tobbb 

25.202/78 

236 

1328 

1. 

7 2 • 4666bb 

25.202/78 

248 

1256 

5 0 " 

72.425000 

25.200000 

2 b 1 

1197 

0. 

72.4«a«ay 

25.183333 

III 

1296 

u 

72.145111 

25.233333 

l IB 

764 

a. 

/2,42777b 

2 5 , 1 8 3 3 3 3 

285 

1206 

9. 

72.175000 

25,222222 

286 

825 

10. 

72.4777/7 

25. lbbbb / 

299 

1288 

11. 

72.022224 

25,21666 / 

328 

607 ’ 

12, 

/2. 022224 

25.211111 

437 

609 

14, 

/2.14194b 

25. 193056 

338 

/ / 6 

1 A. 

7 2.U44444 

25.205555 

340 

644 

i 5, 

/ 2 . 1 / 7 / / 8 

25,1 / / / 78 
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849 

10. 

7 2.22/7 / 7 

25, 14/222 
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i / . 
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ia . 

72.452/ / / 
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1 1 3b 

19. 
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26,119444 
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754 
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460 

675 
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72.441667 

25,058333 

467 

1283 

26. 
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483 

909 

27, 
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734 

26, 

72,005555 

25.105556 

497 

631 
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1237 
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4 2, 
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7 1 . 8 0000 0 
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212/ 
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